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We propose and analyze Maxwell’s demon based on a single qubit with avoided level crossing. Its 
operation cycle consists of adiabatic drive to the point of minimum energy separation, measurement 
of the qubit state, and conditional feedback. We show that the heat extracted from the bath at 
temperature T can ideally approach the Landauer limit of ksT In 2 per cycle even in the quantum 
regime. Practical demon efficiency is limited by the interplay of Landau-Zener transitions and 
coupling to the bath. We suggest that an experimental demonstration of the demon is fully feasible 
using one of the standard superconducting qubits. 


Controllable small systems, e.g. in form of molecules, 
micro-beads or nanoelectronic circuits m, have re¬ 
cently made it possible to experiment on thermodynam¬ 
ics in the regime where fluctuations play an important 
role. One of the most interesting issues in such stochas¬ 
tic thermodynamics is the role of information. On the 
level of thought experiments this discussion dates back 
all the way to Maxwell’s demon (MD) [7] and then to 
Landauer’s principle 0. Till now, majority of the stud¬ 
ies, both theoretically and in particular experimentally 
SHE], have been focusing on the classical regime; see, 
however, nsnzi and references therein for theories on 
quantum systems. The aim of this Letter is to present a 
simple and experimentally feasible quantum MD that can 
operate at the limit of thermodynamic efficiency. Demon 
operation is based on a qubit, a two-level system (TLS) 
with avoided level crossing between the ground ( g ) and 
the excited (e) state as illustrated in Fig. [l] The TLS 
is coupled weakly to a bath at temperature T. The MD 
moves around the level crossing by tuning the control 
parameter q. 

We first present an ideal operation cycle and results 
for the MD without justifying their feasibility. The cycle 
start at point A in Fig. |I] where the level separation 
A E is maximum denoted by A Ea- If A Ea kgT, the 
system is in the ground state with zero entropy. In the 
first leg, q is moved adiabatically to the point X with 
minimum level spacing A Eq- We assume that the ramps 
are linear in time, i.e. q = constant within each leg. 
Ideally, for AE 0 -C kgT and fully adiabatical ramp, the 
entropy of the TLS increases by ks In 2, corresponding 
to heat extracted from the bath equal to Landauer value 
Ql = ksT In 2. This part of the cycle is equivalent to 
that in a classical TLS - see Ref. HU, although the non¬ 
vanishing AE 0 creates essential differences in the quan¬ 
tum regime. On reaching point X, a measurement of the 
qubit state is performed. If the detector measures the ap¬ 
propriate quantity: the energy or the curvature d 2 E/dq 2 
(which, depending on the physical realization of the sys¬ 
tem can represent, e.g., effective capacitance or induc¬ 
tance in an electric circuit) it does not perturb this state 



FIG. 1. The scheme and operation cycle of the qubit demon. 


and therefore, does not incur any thermodynamic costs. 
Note that, in contrast to a classical MD, the possibility 
to realize such a measurement required for the thermo¬ 
dynamically efficient quantum MD, is not automatic, but 
relies on the feature of adiabatic dynamics to leave the 
qubit density matrix diagonal in the energy basis. Af¬ 
ter the measurement, the detector, ideally, “knows” with 
certainty whether the system is in the ground state or 
in the excited state in each cycle. In the repeated cy¬ 
cles, the two states are realized with equilibrium thermal 
probabilities. This makes it possible to apply conditional 
feedback: if the system is in the ground state, the control 
gate is moved quickly back to the original position. In 
the excited state, a 7r-pulse is applied to the qubit swap¬ 
ping the state e —>• g, the process in which the work done 
by the source is — A E 0 HOT, followed again by a fast 
ramp in the ground state back to the original position A. 
Extracted heat, averaged over several such ideal cycles is 

(0) = -r 1 ln(l + ^)+ 1+ ^ a , (1) 

where /3 = {ksT)~ l . 

Several processes can compromise the ideal perfor¬ 
mance of the qubit demon. First, one may extract less 
heat from the bath if the ramp A —i X is not fully adi¬ 
abatic. Also, the feedback ramp X —i A needs to be 
“fast” to avoid bath-induced transitions but “slow” to 
avoid Landau-Zener excitation out of the ground state. 
Both processes would create extra dissipation. To con- 
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sider them quantitatively, we start with the operator of 
the power injected into the system, P(t) = Hit), where 
H(t) is the system Hamiltonian. Work W done on the 
system over a time interval r is then W = fj P(t)dt. 
Since we are mainly interested in quantities averaged 
over many cycles, it is convenient to describe the dy¬ 
namics with the density matrix pit). For average work 
(W) = (JJ P(t)dt) one finds by elementary analysis [20] 

{W) = (H(t)) - (H( 0)) - f Tr [mm] dt. (2) 

Jo 

Equation ([2]) represents the first law of thermody¬ 
namics, combining the change of the internal energy 
and the average heat generated in the bath: (Q) = 
— Tr[p{t)H(t)\dt. In the basis |e)} of instanta¬ 
neous eigenstates of the qubit, this expression gives 



FIG. 2. (Q) as a function of the rate v = hq/AE a in a quasi¬ 
static A —> X ramp. The horizontal arrow points the result 
of fully adiabatic evolution based on Eq. 0- The solid black 
line is the full numerical solution of the master equation Q. 
The blue dots are based on the stochastic quantum trajectory 
simulations. The parameters are e = 0.1, /3AEa = 10, and 
g 2 = HT 3 . 


(Q) = [ dtp gg (t)AE(t). (3) 

Jo 


The standard master equation for the qubit density 
matrix p(t), ignoring pure dephasing, is 

p gg = -melpgePfo AJ5(0*7»] _ Tx Pgg + T l 
p ge = n(p gg - 1/2 )e-*/o*-**<*')*'/» _ ir s p 9e , (4) 

where Q = 3 ^_ e2 q depends on the ramp rate of q , 

e = AE 0 /i2AEA), and the level spacing is A E(t) = 
2AEA^/qit) 2 + e 2 . Point A corresponds to q = —1/2, 
whereas point X is at q = 0. The relaxation and excita¬ 
tion rates of the qubit are, respectively, 


W) 


g 2 a m r 

l_ e -/3A E(t)> ^ 


g 2 A E(f) f 5) 
h - 1 ’ y ’ 


and Ts = F| + F-)-. Parameter g 2 denotes the dimension¬ 
less strength of coupling to the bath. 

To learn about stochastics beyond the averages above, 
we analyze quantum trajectories corresponding to the 
master equation (0) for the evolution of the wave func¬ 
tion | ip(t)) = a(t)\g) + 6(t)|e). According to the quantum 
jump approach [21j . the probability to relax to | g) (get 
excited to |e)) over a short time interval At is given by 
Api = r^|6(t)| 2 At: (A= P-|-|a(t)| 2 At). If such a jump 
does not occur, the system evolves as 

a = AE{tt)df/h b+ ^ATa\b\ 2 

b = +ln e +i ti AE ^) dt '/ n a - ^Ar6|a| 2 , (6) 

where AT(t) = T±(t) — r-j-(i) = g 2 AE(t)/H. The evo¬ 
lution of the system is obtained numerically using a 
stochastic simulation based on these equations. The heat 
dissipated into the bath is given by (—)A E when the sys¬ 
tem relaxes (gets excited) [22] , see also [25] . 


We analyze the MD cycle and errors quantitatively 
assuming that the measurement and the conditional tt- 
pulse are fast enough (duration -C T^ 1 ) so that no heat is 
exchanged in this time interval. Consider first the quasi¬ 
static sweep A —>• X in time r = (2q)~ 1 . For r —> oo, the 
population obeys 

Pgg(t) = Fi/r E = (1 + e-^W)- 1 ( 7 ) 


according to Eq. 0. Inserting this expression into 
Eq. (§ , we obtain 


(Q) 


1 + e 0A E 0 AE 0 

V1 + ed AE ^ J 1 + e~P AE ° 


AE a 

1 + e~d AEA 


( 8 ) 

Assuming everywhere below that (3AEa ^ 1, we recover 
Eq. Q. If AEq = 0, we obtain then the most familiar 
result (Q) = —ksT In 2 = — Ql- For such a degenerate 
case, the average heat equals the heat obtained in each 
adiabatic realization, but this is not so for A E 0 > 0. 

For finite values of r, the evolution becomes non- 
adiabatic and can be analyzed using Eq. 0. Provided 
h/i2ne 2 AEat) <C 1, one can neglect Zener tunneling 
processes and set = 0 in Eq. 0. Then one finds ex¬ 
actly 



dt ,pY Y ,(t)AE(t)AE(t')e-fi dsT ^ 

4cosh 2 [/3AE(t')/2] 

(9) 


In Fig. 0 we present results on (Q) normalized by fcgT 
as a function of the ramping rate v = H/^tAEa) in the 
A —> X ramp for chosen parameters. The black line is 
based on the master equation 0 together with ( 0 , and 
is indistinguishable from Eq. 0. The blue symbols are 
from the stochastic simulation with 50000 repetitions for 
each point. The horizontal arrow indicates the adiabatic 
result of Eq. 0 . 

The errors in the fast return ramp X —> A in time t# 
are produced by the Landau-Zener processes which gen¬ 
erate finite population of the excited state. If the excited 
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state amplitude is bf at the end of the sweep, there is 
“quantum” heat ~ \bf\ 2 AEA dissipated into the bath in 
the next adiabatic ramp to q = 0. On the other hand, 
if the sweep is too slow and T^tr <C 1 does not hold, 
the system can be excited by the bath, and there is ex¬ 
tra “classical” generated heat when the system relaxes 
either before or during the adiabatic ramp in the next 
cycle. These two errors constrict the ramp time to the 
interval h/(e 2 AEa) <C tr -C IT 1 as a condition of ap¬ 
proaching the ideal operation. The lower bound is the 
same as for the standard Landau-Zener dynamics. 

We estimate these two contributions to heat as fol¬ 
lows. For weak dissipation and system being mostly in 
the ground state, we write the evolution of b as (see <§) 
b ss ±fl e + l f 0 AE ( t ) dt / h . Linearizing the expression in 
the exponent with respect to q, and denoting the dimen¬ 
sionless rate as vr = H/(2trAEa), we obtain the final 
population at the end of the ramp as 

\bf\ 2 *v%/{ 16e 4 ), (10) 

valid for \bf\ 2 -C 1. This result agrees with the one ob¬ 
tained in [23] . 

The corresponding classical error is determined by the 
probability p-j- of a thermal excitation g —► e during the 
ramp, again assuming that |n| 2 ss 1. This yields 


Pt 




^g 2 e 3/ 2 e-^ 
VPAEa 


J R ■ 


( 11 ) 


Both excitation processes, described by Eqs. (10) and 


(11), combine to produce extra average heat dissipated in 
a cycle: (Q d iss) ~ ( \bf \ 2 + Pf)AE A . Minimum of (Q diss ), 


with respect to vr is then obtained at 

_ , 80rgW 2 e- 2 ^^ u/3 
VR,opt ^ VPAE a ^ 


(13) 


In this optimization, we considered only the most 
straightforward linear return ramp, for which the exci¬ 


tation amplitude bf (101 decreases as 1 /tr. In principle, 
bf can be suppressed more strongly if the return ramp is 
switched not abruptly, as in the cycles considered above, 
but smoothly. As discussed in the Supplementary Mate¬ 
rial [25], the amplitude bf decreases as (h/AE qtr) tl , if 
only the nth-order derivative of the switching function 
is discontinuous, while all previous derivatives are con¬ 
tinuous, and exponentially, if the switching function has 
continuous derivatives of any order. 

The approximations above are shown together with full 
master equation solutions in Fig. [3] again for e = 0.1 
and PAE a = 10 for the ramp X -> A as a function of 
the sweep rate vr for two cases of g 2 = 10 -3 and g 2 = 0. 
The predictions based on the analytical estimates, i.e. 



FIG. 3. Population of the excited state at the end of the 
return sweep X —> A based on the master equation (U| (sym¬ 
bols), together with the analytical approximation of \bf | 2 + p d 
(lines) for the parameter values e = 0.1, /3AEa = 10, 
g 2 = 10 -3 (upper line and symbols) and g 2 = 0 (lower line 
and symbols). 


1 — p gg j = \bf\ 2 + pf are shown by black and blue lines 
for the two values of g 2 , respectively. The correspond¬ 
ing solutions based on the master equation are shown by 
solid symbols for the two values of g 2 as well. We see 
that the approximations are relatively accurate in this 
regime. The decline for vr —> 0 of the master equation 
data set for g 2 = 10“ 3 (black dots) is due to relaxation 
back from the excited state, the effect ignored in the an¬ 
alytical estimates. 

As the last step of the numerical analysis of the demon, 
we have performed simulations based on the quantum 
trajectory method in cyclic operation. In this situation, 
once a cycle has finished, the state of the qubit when 
reaching point A serves as a new initial state of the next 
adiabatic ramp A —> X. Thus, dynamics of the previous 
return path makes the initial state not necessarily the 
ground state, which leads to extra dissipation. Figure [4] 
shows results of this analysis, where heat accumulation 
over N cycles up to N = 10 4 has been analyzed under 
different driving parameters. In all these simulations, e = 
0.1, PAEa = 10, and g 2 = 10 -3 . In Fig. |4](a), the rate of 
return sweep remains constant, at vr = 1.25-10 -3 , which 
is close to its optimum value. The parameter indicated 
within the figure is the ramp rate v of the adiabatic leg. 
We see clear cooling by the MD and correspondence with 
Fig. For comparison, the fast rising (red) line indicates 
results of the calculation where the parameters are as 
for the lowest black curve in the figure, but without the 
measurement and 7r-pulse before the ramp back to A. 
This protocol is similar to a periodically driven two level 
system, which has been studied in detail recently BSJH7]. 
Naturally in such deterministic cycles there is positive 
heat dissipated into the bath according to the second law. 
The lowest (blue) line shows the result of Eq. |l]) which 
would be achieved if the first ramp were fully adiabatic, 
and the system would remain constantly in the ground 
state during the return leg. In Fig. [4] (b) we show similar 
data as in (a), but now keeping the adiabatic ramp at a 
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FIG. 4. Results of stochastic quantum trajectory simulations 
over N cycles of demon operation. In (a) the quasistatic ramp 
rate v (with vr = 1.25 • 10 -3 ) and in (b) the return ramp 
rate vr (with v = 1 ■ 10 -< ’) are varied as indicated by the 
numbers with arrows. The top (red) line in both (a) and (b) 
shows the result when the demon is not active, whereas the 
bottom (blue) line shows the result of elusive ideal cycles. The 
horizontal dashed line shows the break-even result. 

constant rate v = 10 -6 . The return rate vr is in turn 
varied here. We see correspondence with Fig. [3j showing 
the optimum rate of about vr = 1.25 • 10~ 3 . 

It is possible to make some general statements as to 
the thermodynamic efficiency of the considered quan¬ 
tum MD. Thermodynamic cost of information acquired 
by the detector, based on Landauer’s principle, is 
-k B T(pgg\npgg + p ee In p ee ), where p ee = 1 - p gg . For 
ideal cycles considered above, this cost equals the heat 
|(Q)| ofEq. 0 extracted in a cycle. In this respect, ide¬ 
ally, the presented MD is fully reversible. Work W„ done 
in applying the 7r-pulse can obtain two values, W v = 0 
if the system is in the ground state and nothing is done, 
or W n = —AEo if the two-level system is in the excited 
state and the pulse is indeed applied. Therefore, averag¬ 
ing over many cycles, we obtain 

A P 

{Wn) = -^e A, E 0 = ~ 1 + e p AEo - ( 14 ) 

Inserting Eq. Q for average heat, we find that 

{W w )-(Q)=p~ 1 \n{l + e~P AEo )>0. (15) 


Both (Wtt) and (Q) are negative. Thus the work extrac¬ 
tion in this step is less than the heat extraction in the 
whole cycle. In steady state, the average work in the 
whole cycle (W) = (Q) according to Eq. <1 , see also 

m- 

Finally, a few words about experimental feasibility of 
the proposed qubit demon. The principle studied here 
applies naturally to a large variety of systems, includ¬ 
ing solid-state qubits. As a concrete example, we may 
consider a superconducting qubit, e.g., a charge qubit 
[28l [29] or transmon qubit [291 SO] ■ Typical energy scale 
in this case is AE^/k B ~ 1 K, as aluminium, the most 
common superconductor in this respect, has an energy 
gap A/ks ~ 2 K. Since the operating (bath) temper¬ 
ature of these devices is T ~ 100 mK, we indeed have 
/ 3AEa ~ 10, as in the numerical examples above. In a 
charge-type qubit, e is the ratio of the Josephson cou¬ 
pling and capacitive charging energies, and can be made 
to be ~ 0.1, based on design and tuning by magnetic flux. 
The value of g 2 = 10^ 3 would correspond to a relaxation 
time of the qubit, Tj = 1/Ts, which is then in the range 
0.1...1 /us, quite feasible for the current technology. The 
ramp times for the demon fall also in a favourable regime. 
Since r ~ 10 5 H/AEa ~ 1 gs yields still about 80% effi¬ 
ciency in the adiabatic leg as compared to the infinitely 
slow ramp (see Fig. 0, we may take this as a naive trade¬ 
off value in terms of efficiency vs. power extraction. The 
duration of the whole cycle is determined almost entirely 
by t. It is interesting to see that 1 p s is about 100 times 
faster than a typical electron-phonon relaxation time at 
these low temperatures [3T1 32] . Therefore, if the bath is 
composed of, e.g., a nricron-scale metallic calorimeter, it 
is reasonable to expect that the steady-state temperature 
of this absorber would decrease once the MD is set into 
cyclic operation [35]. Measuring this temperature drop 
would serve as a smoking gun of the proof of the principle 
of the demon. Last, we found that tr ~ 10 3 h/ AEa ~ 10 
ns would be optimal with these parameters of the MD. 
This again is fully compatible with the time scales of the 
gate and flux operations of superconducting qubits. 
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APPENDIX I: SUPPRESSION OF THE 
LANDAU-ZENER PROCESSES BY SMOOTH 
DRIVE 


In this Section of the Supplementary Material, we pro¬ 
vide a brief discussion of the effect of smooth drive on 
the amplitude bf of the Landau-Zener excitation of the 
qubit out of the ground state in the process of rapid ramp 
back from the level-crossing point to the initial state. The 
most basic ramp-back process studied in the main text 
is characterized by the time dependence of the control 
parameter q(t), which is continuous itself, but has dis¬ 
continuous first derivative: q{t) = 0 for t < 0 (we take 
the instant of measurement and, when needed, the re¬ 
pulse, to be t = 0) and q(t) = qt = t/(2Tn) for t > 0. 
This discontinuity of the first derivative leads to the ex¬ 
citation amplitude bf which decreases as 1 /tr with the 
ramp time - see Eq. (10) of the main text. The point 
of our discussion here is to demonstrate that the ampli¬ 
tude bf can be made to decrease faster with the ramp 
time if any discontinuities of the time dependence 
q(t) of the control parameter at the point t = 0 where 
the ramp is switched on, are suppressed more strongly. 
The statement that we want to make here, demonstrated 
by examples, is that the amplitude bf can be made to de¬ 
creases as 1/rjj, if the switching function q{t) is such that 
only its nth-order derivative is discontinuous, while all 
previous derivatives are continuous, and exponentially, if 
q{t) has continuous derivatives of any order. 

To do this, we need to calculate the excitation prob¬ 
ability | bf | 2 for evolution starting at the level-crossing 
point. For the part of the MD cycle, ramp back to the 
initial state, which is considered in this discussion, we 
are interested in the regime of sufficiently slow time evo¬ 
lution, when the probability |&/| 2 is small, i.e., the qubit 
stays predominantly in the ground state. On the other 
hand, this ramp of the control parameter q(t ) should be 
fast enough for the thermal-bath-induced transitions to 
be negligible in the calculation of bf. As follows from 
Eq. (6) of the main text, in this regime, the evolution 
equation of the excitation amplitude b is 

b= ^n e +i ^ AE( - t ’' )dt ' /n , ( 16 ) 

and can be integrated to give the following expression for 
the amplitude at the end of the ramp: 


bf = 



(e/2)dq 2 iAEa 

2 , 2 ex P{- 7— 

q z + e z n 



(17) 

As discussed above, we are interested in the time evo¬ 
lution which is adiabatic with respect to the Landau- 
Zener transitions. Quantitatively, this means that the 
parameters in Eq. should satisfy the conditions 
1 e H/(AE 0 tr). In this case, the first of these 
inequalities ensures that the integral converges at 
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q -C 1, and can be rewritten in terms of x = q/t as 
follows: 



dx 

x 2 + 1 


exp{ 


iAEo 

h 


/**<*>+ !,■'»>. (18) 

Jo V 


Note, however, that by extending the limit of integration 
to infinity, we are neglecting possible contribution to bf 
from the discontinuity in the control parameter drive at 
the switching-off point t = tr. In principle, the time 
dependence of q at that point also needs to be smoothed 
out in the same way as at the switching-on point t = 0 
we are discussing. 

Equation (18) is valid for any time dependence of the 
control parameter q(t) leading from q = 0 at t = 0 to 
q = 1/2 at t = tr. To obtain specific results, we first 
consider q(t) which satisfies this requirement, and at the 
same time, has a discontinuous nth derivative at t = 0: 


q(t) = t n /(T% + t n ), t> 0. (19) 


Since only the small-g part of the time dependence 
q(t ) affects the integral (18), i.e., the part at t <C tr, 
Eq. (23) can be inserted into Eq. (18) in the simplified 
form, q(t) = (t/r^)™. Then y(t ) = nt 71 ^ 1 /(6 Tr) and 
y(y ) = ny 1 “ 1 /"/(e 1/, ™T^). With this expression, the in¬ 
tegral (fl8|) takes the form: 


bf 


1 f°° dx r iAE 0 T R e x ' n 

2 J a iMu exp{ - ids - 

x f %!T |I - 1, “V + l) 1/2 }, 

Jo 


( 20 ) 


and can be transformed further by the substitution u = 
(q/)V n > v = (ex) 1 / 71 : 


, n f°° dvv n 1 r xAEqTr 

b ' = T t l wm+1 exp{ ^r- 

x [ V du[{u 2 7e 2 ) + l] 1 / 2 }. (21) 

Jo 

The integration region in this equation (real half-axis of 
v ) can be rotated into the imaginary half-axis, demon¬ 
strating explicitly that in the adiabatic limit, when 
e 1 /" H/ AEqTr (for n = 1, this is the same condition 
that ensures that the usual Landau-Zener probability is 
small), the integral is converging at v <C e 1 /", and can 
be simplified and evaluated accordingly: 



FIG. 5. Results on Q (black line) and W (red line) of stochas¬ 
tic quantum trajectory simulations over N cycles of demon 
operation. The quasistatic ramp has the rate v = 1 • 10” 6 , 
and the return ramp rate vr is varied as indicated by the 
numbers with arrows. The bottom (blue) line shows the re¬ 
sult of elusive ideal cycles. The horizontal dashed line shows 
the break-even result. The other parameters in the simulation 
are e = 0.1, /3AEa = 10, and g 2 = 10 -3 . Since temperature 
is rather low, the system happens to start in the ground state 
in all the three presented examples, and therefore W either 
coincides with Q or exceeds it by A Ea everywhere. 


This implies that one can make the the amplitude bf to 
decrease exponentially with tr , if the switching-on pro¬ 
cess q(t) has continuous derivatives of any order. Since 
q(t) = 0 for t < 0, this can happen only if q(t) is non- 
analytic at t = 0. As an example, we take 


q(t) = 


— ^A-TR/t 


t> 0. 


(23) 


Estimating the integral of the evolution equation (16) for 


this q(t) with the exponential accuracy by the steepest- 
decent method, we find 


\b f \ ~ e - ( - 2AEoTR/R)1/2 . (24) 


This example shows explicitly that the excitation am¬ 
plitude bf can be made exponentially small for smooth 
switching-on process of the control parameter ramp. 


APPENDIX II: WORK AND HEAT IN CYCLIC 
OPERATION OF THE DEVICE 


<22) 


We see that, indeed, the discontinuity of the nth order 
derivative of the time dependence of the control param¬ 
eter q{t) at the switching-on point t = 0 implies that the 
amplitude bf of the qubit transition into the excited state 
in the ramp-back part of the MD cycle is suppressed with 
the increasing ramp time tr as 1 /tr. 


To substantiate our general claim that in cyclic op¬ 
eration work and heat per cycle are equal on average, 
(IE) = ( Q ), we have analyzed in Fig [5] numerically the 
cumulative work and heat by the same method as pre¬ 
sented in Fig. 4 of the main text, but now only over 
N = 100 cycles. The difference between the two quanti¬ 
ties is the energy stored in the qubit as compared to its 
initial value when the operation starts. In the long time 
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limit (N —> oo), the averages of the two quantities per 
cycle coincide. 

One should bear in mind that the data in Fig. [5]cannot 
be interpreted as a result of many measurements of work 
in a single trace shown. It rather presents the evolution 


of Q extracted from the bath. Then, assuming this heat 
trajectory, work is measured after the iV:th cycle. So 
the work values represent a likely outcome of the mea¬ 
surement along the trajectory, if W was measured at this 
point to terminate the protocol at the given value of N. 


